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1. Introduction 

High-order perturbative calculations of non-relativistic heavy quark-antiquark systems are re- 
quired for precise quark mass determinations from QCD sum rules or the lowest upsilon states 
(bottom quark mass) or the energy dependence of the threshold top-quark pair production cross 
section in e^e^ collisions (top quark mass). Concerning the tt quark cross section the status is 
as follows: the NNLO calculations performed about ten years ago revealed a large uncertainty, up 
to ±25%, in the cross section in the resonance peak region |jl|, ^. Subsequent calculations that 
include a summation of logarithms of a, find a much reduced scale dependence around ±(3 — 6)% 
i, |]- The main effect comes from the logarithms at the third order (NNNLO) rather than the 
all-order series [^. Since at NNNLO the ultrasoft scale nitaj ~ 2GeV appears for the first time, 
a complete calculation of the (non-logarithmic) NNNLO correction is needed. Similar conclusions 
apply to bottomonium systems, with larger uncertainties. 

Recently the NNNLO correction to the S'-quarkonium wave-functions at the origin (corre- 
sponding to the residues of the poles of the heavy quark current spectral functions) from potential 
insertions and ultrasoft gluons have been completed [|6[ 0]. The talk presented at this conference 
summarized these results together with NNNLO results on the full energy-dependent spectral func- 
tion ^ relevant to the tt cross section. Since the combined result of [^, ^] has already been dis- 
cussed in another proceedings article we focus here on the full spectral function and the case 
of the top quark. That is, we consider the two-point function 

{q''q'-g'''q^)n{q^) = iJd'xe'^^a\Tif{x)f{Om) (1.1) 

of the electromagnetic top-quark current = ty^t, choosing q^ = (2m, ±£,0) with the pole 
mass of the top quark and E of order of a few GeV. The width of the top quark is taken into 



account by simply letting E E + iVt become complex []10|]. However, one should note that 
a fully consistent treatment of the top-quark decay beyond the NLO approximation requires the 
inclusion of many other electroweak effects that ai^e not yet known. 



2. Remarks on the calculation 

After integrating out the hard and soft momentum scales the problem is reduced to the calcu- 
lation of a non-relativistic correlation function G{E) to third order in non-relativistic perturbation 
theory. The perturbations consist of potential insertions (instantaneous, spatially non-local oper- 
ators) and ultrasoft gluon interactions with the top quarks. Since an infinite number of potential 
(Coulomb) gluons can be exchanged between the heavy quarks without parametric suppression, 
the free heavy quark-anti-quark propagators are promoted to the Green function of the Schrodinger 
operator Hq = —V'^/mt — {asCf)/r with the colour Coulomb potential (Cf = 4/3). Computing 
Feynman integrals with Coulomb Green functions while simultaneously regulating all divergences 
dimensionally to be consistent with fixed-order matching calculations is the main challenge of the 
NNNLO calculation. The third-order correction to G(£') is composed of 

53G = -{0\GodVyGodViGodViGo\0)+2{0\GQdViGodV2Go\0)-{0\GodV3Go\0) + dGus, (2.1) 
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where Go = (^o — E — ie)^^ , |0) denotes a relative position eigenstate with eigenvalue r = 0, and 
5Vi the ith order perturbation potentials (see [§]). The contributions to G{E) involving only higher- 
order corrections to the Coulomb potential have already been computed in [ jlTl ] and are included in 
the following numerical result. The ultrasoft contribution is (D = J — 1) 



8Gu 



-i){igsfCF 



X 1 



X 1 



—{in) 5( \p,-P2) + {ig.) Y^^^r^ 



nit 



m 2 (PA-ps)^ 



iG^^\p3,PA;E + k°) 

iG^'\p5,P6;E) (2.2) 



with &^-^\p,p';E) the colour-singlet/octet momentum-space Coulomb Green functions. In posi- 
tion space this expression simplifies to three instead of seven loop integrations and similar simpli- 



fications apply to the other terms in ( |2.1| ). However, the integrals are divergent and the 1 /e poles 
must be extracted in momentum space. A further complication is that the Coulomb Green functions 
are not known in D dimensions. The strategy therefore consists of identifying all divergent sub- 
graphs, and to calculate them in (i-dimensional momentum space. Then combine the result with the 
sub-divergence counterterms related to the renormaUzation of potentials and non-relativistic cur- 
rents and perform the remaining integrations in three dimensions. Due to the incomplete treatment 
of finite width effects an over-all divergence a/'^/e x F, remains in the tt cross section, which is 
minimally subtracted in the result below. 



3. Size of logarithmic and non-logarithmic terms 

We consider the residue of the correlation function at the lowest-energy bound-state pole at 
El = — wjf (ttvCf )^/4 + ... to compare the new NNNLO non-logarithmic terms [Q, ^ to the previ- 
ously known logarithms [12, |l^, 14, 15], since in this case a simple numerical result can be given. 
Zi , defined by 

2n e^e, 3 Zi 



(3.1) 



2mJ Ei-E-ie' 

is related to the height of the cross section peak by Rpcak ~ 187re^Zi/(/M,^r,), so we expect similar 
conclusions to hold for the entire tt cross section. The NNNLO expression for Zi reads 



871 



1 + a, -2.13 + 3.66L 
+ af 8.38-7.26 Ina,- 13.40L + 8.93L- 
+ 11.01 + [31.5S]c3,nf - 9.79 In a, - 16.35 In^ a, 
+ (53.17-44.27 lna,)L-48.18L2 + 18.17L^ 

{mtasCp)^ 



871 



x ( l-2.13a, + 22.64a^2 + [-32.96 + [37.58]c3,„^]a,^ ). (3.2) 
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Since the scales m^, mja^ and mta^ are all relevant here, there is an ambiguity in the representation 
of the logarithms. The above result uses = o;s (/x) and puts the explicit /i -dependence into L = 
In jj. /{mtCFOCs). The NNNLO result is not yet complete: the constant 1 1.01 includes an estimate of 
the NNNLO connection to the Coulomb potential, 03 = 3840 [16], and sets certain unknown ^{e) 
potential terms to zero. This is expected to have a minor effect [^]. More important is that only the 
«y^-parts of the third-order matching coefficient C3 of the non-relativistic current Y^^^'X known 
[17], which turn out to be very large ([37.58]c3,,iy). In the following numerical results for the cross 
section we therefore consider two options, one where the constant part of C3 is set to the known n f 
terms, the other where it is set to zero. (The logarithms are all known and always included.) 

We observe (third and fourth line of (3.2)) that the typical size of non-logarithmic terms of 
individual third-order corrections (ultrasoft, non-Coulomb potentials, Wilson coefficient) is about 
40a^"' 10% (> 100%) for toponium (bottomonium). However, large cancellations between indi- 
vidual terms and between logarithmic and non-logarithmic terms occur. Thus, to obtain a reliable 
third-order result the non-logarithmic terms are crucial and a final assessment needs the missing 
nj-independent term in C3. This is seen in the last line of (3^), which shows the series for a, = 0.14 
where L = 0. The large NNLO correction is evident. On the other hand, the third-order correc- 
tion is not anomalously large, although the final coefficient will only be known when the term 
[37.58]c3,«^ is replaced by the full result for C3. The NNNLO result shows a strong reduction of the 
scale dependence compared to NNLO as discussed in [^, but the perturbative prediction becomes 
unstable for /i < 20GeV. A study of this problem for the Coulomb corrections, where a resum- 
mation of the perturbative series can be done by means of a numerical solution, has shown [11] 
that the perturbative prediction for > 25 GeV is close to the true result, hence we do not consider 
scales n < 25 GeV. 



4. Top- quark cross section 

We next discuss the tt production cross section in e^e^ annihilation near threshold. More 
precisely, we consider the /?-ratio 

R = Onx/oo = l2ne^\mU{q^) {oo = ^naU {3s)) , (4.1) 

neglecting the axial- vector contribution from Z-exchange for the purpose of discussing the impact 
of the QCD NNNLO correction. The top quark pole mass should be avoided as an input parameter. 
Here we use the potential-subtracted mass [ |T8[ ] implemented as explained in [pi]]. The parameters 
for the cross section calculation are: m, ps(20GeV) = 175 GeV, F, = 1.4GeV, a,{Mz) = 0.1189. 

The successive LO ... NNNLO approximations to R are shown in Figure |] (top panel). At 
/I = 30 GeV the size of the third-order correction is up to 10% depending on the assumption for 
C3. When all known terms in C3 are included the peak cross section is about 10% larger than in 
the renormalization-group improved NNLO calculations [^, ^, ^ due to the sizeable constant term 
related to 11.01 + [37.5 8] C3,„^ in (^^. Contrary to the NNLO approximation the third-order result 
shows good convergence of the perturbative expansion. The bottom panel of Figure ^ consequently 
displays a strong reduction of the scale dependence from NNLO to NNNLO. The residual scale 
dependence is similar at NNNLO and in the renormalization-group improved calculations, which 
already captures correctly the logarithms of ^u. It therefore appears that with a complete NNNLO 
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Figure 1: (top) Successive approximations to the 7?-ratio at fixed ji = 30GeV. At NNNLO two implemen- 
tations of C3 are shown as discussed in the text, (bottom) Renormalization scale dependence at NNLO and 
NNNLO. Here all known terms in C3 are included. 



result and a sununation of higher-order logarithms at hand, the demands on an accurate theoretical 
prediction of the cross section near threshold can be met as far as QCD corrections are concerned. 
In particular, the scale dependence of the peak position which is indicative of the accuracy of the 
top mass measurement is now well below 100 MeV. 



5. Summary 

The NNNLO QCD correction to the tt cross section near threshold is now nearly complete. 
We presented for the first time the result of the third-order potential and ultrasoft correction to the 
non-relativistic heavy-quark correlation function. We find that the third-order correction behaves 
well (contrary to the anomalously large effect at NNLO) and removes a large part of the theoretical 
uncertainty. The new non-logarithmic terms are numerically important and increase the cross sec- 
tion relative to the renormalization-group improved NNLO result by about 10%, when all presently 
known terms of the three-loop matching coefficient C3 are included. Further work is necessary on 



a consistent treatment of electroweak and finite-width effects (see [19, 20, 21]). 
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